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Abstract 

One-dimensional motion of a quantum point particle is usually described by 
its wave function ip(x) where the argument x G M represents a (measurable) 
coordinate and where the integrated probability density is normalized to one, 
J if)*(x)ijj(x) = 1. The direct observability of x may be lost in VT— symmetric 
quantum mechanics where a "smeared" metric kernel Q( x ,xf) 7^ ${ x ~ x> ) 
may enter the double-integral normalization, J J ip*{x)Q( x ^ x >^(x') = 1. We 
argue that such a formalism proves particularly suitable for the introduction 
of a nonvanishing fundamental length 9 > which would characterize the 
"smearing width" of the kernel Q( x ,x')- The technical feasibility of such a 
project is illustrated via a toy family of Hamiltonians H^ N '(X) taken from 
paper I (M. Znojil, Phys. Rev. D 78 (2008) 025026). For each element 
of this family the complete set of all the eligible metric kernels 0/ n(A) is 
constructed in closed form. We show that at any pre-selected non-negative 
fundamental length these metrics can be made to vanish unless \x — x'\ < 9. 
The strictly local inner product of paper I recurs at 9 = while the popular 
CVT— symmetric option requires 9 = oo in this language. 



1 Introduction 



The introduction of a minimal length scale (denoted, say, by symbol 6) is 
tempting on empirical as well as purely pragmatic grounds. Typically, its 
existence would facilitate the regulation of the high energy asymptotics in 
field theory p], etc. A more ambitious motivation of its introduction might 
be sought in M theory or string theory in their various limiting cases [2]. 
Connections between 6 > and the emergence of certain singularities with 
nontrivial physical meaning could further be sought in cosmological applica- 
tions of quantum theory [3]. In the astrophysical context, last but not least, 
the fundamental length might be identified as accounting for the dark energy 
[I] or for the inflationary era in the early evolution of the universe [5J. 

In the majority of similar considerations the fundamental length scale 
emerges as a free parameter. For the particular quantum dynamics of space- 
time, for example, its value can be related not only to the Planck length but 
also, say, to a nonvanishing cosmological constant or to the vacuum energy 
density [6] . In the simplified context of quantum mechanics this quantity can 
even be treated as one of phenomenological, experimentally determined char- 
acteristics, say, of a condensed- matter system [7]. In parallel, in some more 
ambitious theoretical studies the introduction of a fundamental length con- 
stant is being based on a deeper principle like the stabilization requirement 
imposed upon relativistic algebra [8]. 

Many of the latter ideas are implemented using the assumption of non- 
commuta-tivity of coordinates (cf. a rather nonsystematic selection [HJ [9] of 
some sample references). In what follows we intend to develop a different 
theoretical concept in which the existence of fundamental length will find 
its origin and realization via analytic rather than algebraic considerations. 
We shall particularly be guided by the recent innovative interpretations of 
certain analytic potentials characterized by their so called VT— symmetry 
(for a review, paper [10J is recommended). 

More specifically, in our preparatory section [2] we shall define V and T via 
a sample Hamiltonian H (paragraph 12. ip and explain why we believe that in 
quantum theory of similar models one of the most natural definitions of the 
length-scale 9 > could be based on a suitable particular realization of the 
physical Hilbert space 7i. We show, in paragraph 12.21 how a "smearing of 
coordinates" emerges as allowed by the well known ambiguity of the physical 
inner products in TC. This enables us to conjecture that in general, one should 
be able to remove or at least suppress this ambiguity by the requirement that 
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the range of the smearing of coordinates acquires precisely the pre-selected 
non-negative value 9. 

For a quantitative understanding of the similar quantum models with 
built-in scale 9 a family of very specific illustrative examples is introduced in 
section [3j The feasibility of their analysis is achieved not only by the use of a 
non-perturbative technique of solving Schrodinger equations (based on a dis- 
cretization of coordinates, cf. paragraph 13. II) but also by the choice of a very 
elementary, next-to-trivial interaction (cf. ref. [TTJ or eq. (flOj) in paragraph 
13.21 below). In this way all our Hilbert spaces become finite-dimensional and 
the Hamiltonians become represented by certain 2K— dimensional matrices 
H^ 2K >{X) where the real parameter A controls their non-Hermiticity. At the 
minimal Hamiltonian-matrix dimension 2K = N = 2 this renders all the 
eligible "physical" inner products in Ti available in closed form (paragraph 
13.31) making the discussion of the fundamental length 9 trivial. 

The N = 2 conclusions encouraged us to develop and employ a linear- 
algebraic algorithm of the construction of metrics applicable at any even 
Hilbert-space dimension N = 2K. The symbolic-manipulation results are 
sampled, in section HJ at the two subsequent integers K = 2 and 3. The avail- 
ability as well as unexpectedly transparent matrix structure of the resulting 
matrices of the metric proved crucial for our present fundamental-scaling 
purposes. In essence we revealed that the set of 0^ appears composed of 
subsets in which the matrices acquire a band- matrix structure. The el- 
ementary length 9 (i.e., the range of the smearing of coordinates in the inner 
product) is then identified with a measure of non-diagonality (i.e., with the 
ratio between the number of diagonals and dimension) of these matrices. 

In section [5] we turn our attention, first of all, to the next few higher 
dimensions iV > 8 at which the computer-assisted brute-force determination 
of all of the eligible "physical" inner products in 7i still remains feasible. 
Of course, the very enumeration of the N— parametric sets of the resulting 
matrices 0^) becomes clumsy. For this reason we developed a recurrent 
technique of their description which offers their compact classification as 
well as their explicit and compact description at any dimension iV = 2K. 

The phenomenological core of our present message lies in several observa- 
tions formulated in several separate subsections of section [6j We point out, 
i.a., that our present choice of the family of Hamiltonians H^ 2K \X) looks 
particularly suitable for illustrative purposes since the bound-state energies 
remain real in a dimension-independent interval A G (—1, 1). One finds out 
that at a fixed dimension iV = 2K the phenomenologically motivated fun- 
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damental length 6 can equally well be perceived as a means of classification 
of eligible metrics. This feature of our solvable model becomes particularly 
useful when a fixed choice of the fundamental length 6 > is analyzed in the 
continuous-coordinate limit N — > oo. 

Our last, summary section [7] re-emphasizes the importance of the exact, 
non-perturbative solvability of our schematic benchmark example which ad- 
mits the exhaustive and exact construction of the complete menu of matrices 
@( N ) of the metrics at a given N. In the future, having such a transpar- 
ent methodical guide at our disposal we may expect that all the prospective 
transitions to the more realistic interaction models will be facilitated and/or 
more easily mediated via approximate (e.g., perturbation-theory) techniques. 



2 Nonlocal inner products 

2.1 VT— symmetric models 

Our main source of inspiration can be traced back to the discoveries of 
existence of the real spectra of energies generated by non-Hermitian one- 
dimensional Hamiltonians [T2j [T5]. For illustration let us recollect just the 
Buslaev's and Grecchi's (BG, [13]) Hamiltonian H = p 2 + V(x) ^ W of this 
type containing the asymptotically quartic "wrong-sign" potential 

V(x) = V {BG \x) = -(x-isf + 0(x 2 ) V*(x), x 6 (-00,00), e>0. 

The distinguishing feature of this non-Hermitian Hamiltonian (studied, later, 
also by Jones et al [H]) is that it is characterized by its VT— symmetry 
VTH = HVT where the operators V and T stand for the parity and time 
reversal, respectively. 

The essential merit of the BG model is that the asymptotically dominant 
parts of the general solutions ip\ 2 ( x ) °f ^ ne related differential Schrodinger 
equation near the respective endpoint coordinates x\ t 2 = ±00 are easily de- 
duced, 

4f\x) = c^ 2) ^ G \x) + c^ 2) ^ G \x) , ^ G \x) = gii^W+ow 

As long as e > we can set c + = guaranteeing that the resulting functions 
ip[ B 2 G \x) will asymptotically vanish, ip[ BG \+oo) = and ipi BG \— 00) = 0. 
After analytic continuation their matching near the origin gives the physical 
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bound-state solution which is analytic at all x G (—00, 00) and quadratically 
integrable, i.e., ip^ G \x) G Ej 2 (1R) or 

/ dx W G \x)X ^ n BG \x) < OO. (1) 

It has also rigorously been proved in [13] that the energies E = E n are all real, 
non-degenerate and growing with the main quantum number n = 0, 1, . . .. 

From our present point of view, it is more important that the related wave 
functions remain mutually non-orthogonal. This means that the physical 
information carried by these wave functions remains unclear. In order to 
restore the physical probabilistic interpretation of such a VT— symmetric 
model one must modify the inner product [TQl US]- Usually, this goal is 
achieved by the replacement of the unphysical Hilbert space E, 2 (]R) := HS F > 
(where F stands for "first" or "friendly" or "false" ) by its amendment 7^ 
K, 2 (M) where S means "second" or "standard" and where the bound states 
become mutually orthogonal [IB] . 

2.2 Fundamental length 9 as a measure of non- locality 

During the return to the "standard" Hilbert space one reveals that our 
illustrative potential V^ BG ^ (x) is an extremely "user-friendly" interaction. In 
this sense many of its properties appear rather exceptional (cf. [13] where 
one finds that the transition to Ti^ leaves the interaction local). For this 
reason it makes sense to recall also several other illustrative models. In order 
to make the picture comparatively complete one must recollect, e.g., the 
imaginary cubic oscillator of Bessis and Zinn- Justin [17] and/or the whole 
one-parametric family of its generalizations with V(x)/x 2 = (ix) 5 where 5 > 
(cf. [10] for more details). In all of these cases the transition to makes 
the interaction strongly nonlocal (for illustration we recommend to check the 
details, e.g., via their perturbative illustration at 5 = 1 in paper [18J). 

Let us emphasize that the necessary condition of the possibility of the 
transition from the unphysical Hilbert space TiS F ^ to its physical parallel 
lies in the reality of the spectrum of the "non-Hermitian" Hamiltonian 
in question. Once this condition is satisfied, the most common realization 
of the correspondence TL^ — > usually (cf. [TU]) proceeds under the 

assumption that the space remains spanned by the same wave functions 
and that only the inner products are re-defined as non-local [19]. Here we 
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shall follow the same recipe. Without getting too deeply in the underly- 
ing mathematics let us only recollect that for all the sufficiently elementary 
Hamiltonians the inner product used in TiS s ' may be understood as leading 
to the generalized, double-integral orthonormalization rule 

/ dx dx' ip* m {x) 0(3^/) 4> n {x') = S mn (2) 
Jm 2 

where 5 mn is Kronecker symbol. For the Hamiltonians which are "triv- 
ially" Hermitian in Ti^ the standard textbook scenario characterized by 
the Dirac's metric B( x , x ') — ~ x> ) leads merely to the degenerate version 
eq. ([T]) of eq. fl2J). Vice versa, Dirac-non-Hermitian Hamiltonians H ^ W 
with real spectra necessitate a selection of a Hamiltonian-dependent metric 
kernel Q( x ,x') 7^ 8(x — x') in eq. (j2J). One just replaces the elementary Dirac's 
Hermitian conjugation 

T (F) : ijj{x) -> %l)*(x) 
by its non-local, more complicated version 

rW:V(a:)- / dz,l>*{z)Q M . (3) 
Jm 

There exist examples (based on a sufficiently elementary choice of Hamilto- 
nian H - cf., e.g., [20]) where the metric Q( X}X ') itself can even be constructed 
exactly, non-perturbatively . 

In this language our present key message is that our equations (J2J and (j3J) 
may be complemented by the phenomenologically motivated limited-range 
constraint 

(x , x . o ^ only if \x -x'\<6 (4) 

using any preselected real quantity 9 > 0. One must keep in mind that this 
quantity is just an elementary upper estimate of non-locality imposed upon 
the metric in TiS s \ Its size may be perceived as a measure of the "smearing" 
of the coordinate x which is due to the loss of the direct physical meaning 
and measurability of the real variable x in 

Strictly speaking one should not even call quantity 6 > a "length" . At 
the same time, one feels that once the range of the smearing of the metric 
kernel is assumed restricted by eq. PJ, the effect of this smearing will quickly 
decrease when the measured distances exceed the preselected "fundamental" 
value 9. This observation may be also read as a core of our present project 
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specifying a class of models where the smearing is guaranteed to be safely 
short-ranged in scaled units. 

In this setting the mathematical questions emerge which concern not 
only the existence of similar models (this question will be answered here 
affirmatively and constructively) but also their further properties. In this 
sense, our present paper should be perceived as the first step towards a more 
extensive theory. Certainly, it will be non-Hermitian in (where the 

variable x is only a non-observable auxiliary quantity) but safely Hermitian 
in Ti.^ (of course, only here the fully consistent concept of distance can be 
defined). In this context, equation (j3J) might acquire its proper meaning as 
mediator of coexistence between certain asymptotic locality and short-range 
non-locality of certain less standard models of quantum dynamics. 

3 Toy model 

The choice of the value of parameter 9 must stay compatible with our in- 
tuition and available experimental evidence. For example, even the extreme 
choice of a very large 9 may be tolerated in the context of bound states 
where the observable range of non-locality will effectively be limited by the 
exponential decrease of wave functions [21]. In contrast, for a consistent de- 
scription of scattering one must necessarily require that the bound 9 upon 
the non-locality in eq. (j3j) must at least asymptotically be very small [TT] . In 
this sense it is rather encouraging that there exist several different solvable 
models of scattering where the smearing size 9 in eq. (j3j) is strictly equal 
to zero [Til [22] . The methodical importance of the latter family of illustra- 
tive examples is further underlined by the fact that for the vast majority of 
quantum models with nontrivial metrics the numerical value of 9 happened 
to remain infinite [2TJ [23] . 

On this background one must be careful with predictions of the existence 
of Hamiltonians and metrics for which our upper-estimate 9 of the built-in 
non-locality proves nontrivial, i.e., non- vanishing and not too large. Thus, 
we have to offer an explicit, schematic, constructive example of such a Hamil- 
tonian H and of its specific short-range-smearing metrics 0. 



7 



3.1 Runge-Kutta lattice of coordinates 

Let us replace the differential form of a given Hamiltonian H = p 2 + V(x) 
by its discretized Runge-Kutta approximation leading to the linear difference 
Schrodinger equation 

- " KXM> ~ ^'If + * M + V( Xk ) ii Xt ) =E^„) . (5) 

Using a suitable, finite or infinite cutoff L > we set 

h h 

X-k = ~L>, x_ K+ i = -L + h, ...,x = -~, xi = -, xk+i = L . (6) 

The lattice spacing h = 2L/(2K + 1) decreases with = 2K in both the 
bound-state phenomenological regime (where L should be kept constant, cf. 
paper [21]) and the scattering regime (where L = L(N) should grow with N, 
see paper [TT]). 

On each level of precision 0(h) and for virtually any strictly local po- 
tential V(x) the latter recipe makes the practical numerical solution of non- 
Hermitian Schrodinger equations at a fixed lattice-point-distance h decisively 
facilitated [25]. The discretization of the coordinates reduces also the above- 
mentioned double-integral orthonormalization rule ([2]) to its discrete ana- 
logue so that the integral kernels Q( x , y ) become replaced by matrices 6^2 (A). 
Under suitable mathematical assumptions [19] these metric matrices define 
the inner product between any two elements ip and <fi of our Runge-Kutta 
version of the physical Hilbert space of states H {s \ 

K K 

£ nx n )e^(x)<t>(x nl ):=((m (?) 

n=-K+l n'=-K+l 

(cf. [16] for more details). 



3.2 Minimally non-local interactions of ref. fll 



In our recent studies of scattering [TT], [221 [25] we revealed that the finite- 
range constraint can be satisfied (and that one can even easily reach its 
lower bound 6 = 0) provided that a non-locality is admitted in the potential. 
A "minimal" generalization of this type leads to the following Runge-Kutta 
Schrodinger equation, 

ip(xk+i) - 2ip(x k ) + V>(a*-i) . T/ / v., N 

72 1" V(x k ,X k+1 )lp(x k+ i) 
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+ V(x k , x k ) if)(x k ) + V(x k , x k _i) ^(sa-i) = E ip(x k ) . (8) 

In principle, eq. (jHJ) must be complemented by asymptotic boundary con- 
ditions. Keeping in mind, nevertheless, that the analysis of the scattering 
scenario has already been performed in ref. [TT], we intend to deal with the 
bound-state option only, 



ip(x- 



Kj 







ip(x K 



+1; 



0. 



(9) 



Let us pick up the most elementary nonlocal interaction as recommended 
in ref. [11] and return to the related eq. (jSJ) where just the two coupling 
constants will be different from zero, 



V(x , a?i) = -A , V(x 1 ,x )=+\. 



(10) 



At each finite K = 1, 2, 3, . . . or N = 2, 4, 6, . . . the resulting Schrodinger 
bound- state eigenvalue problem (jEJ) + flSJ) + ( flOl) will degenerate to the diag- 
onalization of the respective finite matrix H^ N \\) with tridiagonal structure, 
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(11) 



(12) 



(13) 



Qualitatively this family of Hamiltonians can be interpreted as a set of dis- 
crete analogues of the exactly solvable VT— symmetric square well with a 
short-range non-Hermiticity [27]. 

The simplicity of our present family of toy Hamiltonians numbered by 
their finite matrix dimensions N = IK = 2,4,... can be perceived as the key 
benefit resulting from our preference of the non-perturbative Runge-Kutta 
discretization method. One can certainly expect that whenever needed, the 
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present 9 > techniques and constructions will remain applicable also to 
some other, less artificial and more phenomenologically oriented interaction 
models. Such a transition to more complicated models has already been 
shown feasible, in [22J, at 9 = and N = L = oo. Similarly, some more- 
parametric models were shown tractable by the same method in ref. [11] 

3.3 Two- parametric family of metrics Q^ N \\) at N = 2 

At N = 2 and A = cosy? closed formulae are available not only for the 
energies E = E± = 2 ± simp but also for the norms of the eigenstates ip±. 
The Hamiltonian if ( 2 ) (cos <£>) nicely illustrates the subtle difference between 
its right eigenvectors \ip±) and their left-eigenstate partners ((^±| at the same 
energy (the latter row vectors are denoted by doubled bras as in [16]), 

~ ( vr; ) • ^'(^d == i« - ( vr; ) • (i4) 

A biorthogonal basis can be formed of these partner eigenvectors. Thus, the 
manifestly non-Hermitian matrix H^ 2 \cos(p) can be reinterpreted as a ma- 
trix which becomes Hermitian in the ad hoc, Hamiltonian-dependent Hilbert 
space of states endowed with a nontrivial Hermitian-conjugation oper- 
ation of eq. ©. 

The key merit of our N = 2 example can be seen in the straightforward 
availability of all of its admissible metrics which vary with two free parame- 
ters t± da nu, 

= 0( 2 )(cos^) = + \1>-))t- • (15) 

The guarantee of the necessary positivity of these metrics reads t± > 
and holds also, in the similar decoupled form, at all the higher dimensions 
N > 2. After the insertion of eigenvectors (fT4|) in (fT5l) we arrive at our first 
fully explicit matrix formula 

q / (1 — cos y?) 2 (t+ + t ) (1 — cosy?) sin<^(— 1 + + t_) 
I (1 — cos</?) sin </?(— 1 + + t-) sin 2 <p(t + + t_) 

Its inspection reveals that up to an irrelevant overall factor it may be re- 
written strictly equivalent superposition 

0( 2 )(A) =a 1 M 1 (2) (A) +a 2 Mf } (A), A = cos<^ (17) 
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with the two new real free parameters a.\ oc t + + £_ and a 2 oc (— 1 + + £_) sin 99 
and with the following pair of manifestly A— dependent sparse-matrix coeffi- 
cients, 



Mf\\) 



1 - A 






1 + A 



M 2 (2) (A) 



1 

1 



'11 



Such a re-parametrization modifies the overall multiplication factor in but 
it still leaves the positivity constraint very transparent, 



Oi\ > 0, 



ai(i - y 



> a 



2 ■ 



N 



(19) 



We merely have to choose any a 2 from interval (— ol\ sin 99, a.\ sinyj). 

The mutual coupling between ai and a 2 is the price to be paid for the 
simplification of the A— dependence of the metric. At N = 2, fortunately, the 
requirement of a band-matrix form of implies that we have to set a 2 = 
so that the positivity of the metric will be guaranteed by the elementary 
inequality a% > 0. At all the higher N = 2K > 2, similarly, the positivity 
of the metric will trivially be guaranteed by the set of requirements a.\ > 
and a 2 = as = . . . = ax = 0. Although the domain of the positivity of 
the matrix 0W(A) will be perceivably larger at iV > 2, its strict boundary 
can only be determined numerically in general (cf. also refs. [2B] for a very 
explicit sample study of the boundaries of the domain of positivity of the 
metric in the space of parameters). 



4 Computer-assisted construction of all N met- 
rics 6^(A) at N = 4 and N = 6 

The mathematical study of the similarity relation 

QH = WQ (20) 

between a Hamiltonian-type operator H and its adjoint H* dates back to 
early sixties [29]. In physics, the first use of such a feature of a sufficiently 
nontrivial and realistic Hamiltonian H ^ H* emerged much later [TJ] . In the 
so called VT— symmetric quantum mechanics [TU| an additional constraint 
has been accepted by which the metric is factorized into a product of parity 
V and the so called charge C or quasiparity Q [TUl [30] . 
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In our present paper we shall simply use eq. (I2"0"j) as an (implicit) defi- 
nition of all the eligible metrics 9 = 0(H). Our computer-assisted method 
of solving this linear set of algebraic equations for the matrix elements of 
G will be straightforward, incorporating also all the standard requirements 
imposed upon the metric and listed, say, in [19J. A priori we shall not as- 
sume the existence of any other observable like charge or quasiparity. Hence, 
in our constructive considerations at finite dimensions only the necessary 
Hermiticity G = and positivity G > of the metric must and will be 
required. 

4.1 Ansatz at N = 4 

Hamiltonian H^(X) of eq. (112p offers the first nontrivial simulation of the 
non-Hermitian dynamics which is purely kinetic near its "distant" boundaries 
±L and which becomes dynamically nontrivial in the vicinity of the origin. 
The coupling A merely connects two points in the middle of the lattice. The 
four eigenvalues of matrix H^ 4 '(X) read 

E± i± = 2 ± i ^6 - 2 A 2 ± 2 V5 - 6 A 2 (21) 

and remain real in the same interval of couplings A G (—1,1) as above. Sym- 
bolic manipulations on the computer enable us to find all the corresponding 
matrices of the metric G^(A), 

" otf (1 — A) ajg(l-A) a 3 a 4 

a s (1 — A) ati (1 - A) + a 3 (1 - A) a 4 + a 2 (1 - A 2 ) a 3 

a 3 a 4 + a 2 (1 - A 2 ) a>i (1 + A) + a 3 (1 + A) a 2 (1 + A) 

a 4 a 3 a 2 (1 + A) a t (1 + A) . 

They may be interpreted as the following sum with four variable real coeffi- 
cients, 

e< 4 >(A) = Gg )a2ia3ia4] (A) = ai Mi + a 2 M 2 + a 3 M 3 + a A M 4 (22) 

where each component is a sparse matrix carrying a specific A— dependence, 

1-A 

1-A 1-A 2 

1-A 2 1 + A 

1+A 



Mi 



1 - A 

1-A 

1 + A 

1+A 



Mo 



12 



Ma 



10 

1-A 1 

1 1+A 
1 



M A 



1 
10 
10 
10 



(23) 



The first three items may also be treated as band matrices, i.e., as a diagonal, 
tridiagonal and pentadiagonal matrix containing merely one, two and three 
nonvanishing diagonals, respectively. 



4.2 Positivity constraint 

As we already indicated in paragraph 13.31 the diagonal metric ®rio 00 i(A) = 
remains safely positive definite inside the open interval of A G 
(— 1, 1), with two plus two doubly degenerate eigenvalues [i\ 2 = 1 — A and 

11^2 = 1 + A. The remaining three matrices in eq. fl23l are indefinite. Each 
of them possesses a pair of positive and a pair of negative eigenvalues which 
are also easily obtainable in closed form. 

Once we decide to fix ol\ — 1 and treat the remaining three parame- 
ters a 2 , a 3 and a 4 as small perturbations, we may also easily establish an 
allowed range of these perturbations for which the positivity of the metrics 
Si 4 ' ,(A) remains robust and guaranteed. 

Of course, starting from N = 4 it is much less easy to describe the 
strict position of the A— dependent boundary dV of the whole (open) domain 
T> of our four real parameters a,\ (= 1),«2,«3 and «4 in which the metric 
fl22|) is positive definite. At this boundary we may expect that the function 
F := det ®\^ a2 a3 a4 i(^) (which is equal to the product of the four eigenvalues 
of the metric in question) will vanish so that our specification of its zeros 
it needed. This task becomes particularly interesting in the maximally non- 
Hermitian dynamical regime, i.e., say, at the couplings A = 1 — e 2 where 
the real variable e remains very small and where the (real) energy levels of 
eq. (I2~T|) get, pairwise, almost degenerate, with £7 ,i = lT £ / 'V% + (e 2 ) while 
E 2 , 3 = 3T£/V2 + 0(e 2 ). 

In the zero-order approximation O (e°) we reveal, by direct computations, 
that the determinant F will vanish whenever a 3 = ±014. For illustrative 
purposes, let us, therefore, accept this restriction to an exceptional subspace 
of parameters and choose the upper sign for the sake of definiteness. The 
same argument applied in the next-order approximation O (e 2 ) leads to the 
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specification of the next quantity a 2 = 1 + and leaves just the single 
parameter in the metric unspecified, a 4 := 7. 



/" Y 


1 \ 






-3- 




x> -5- 





Figure 1: The boundary curve 7 = 7(e). 

Along the boundary dT> we must have F = -F(7, e) = 0. This is an equa- 
tion which establishes an implicit polynomial relationship between 7 and e, 
i.e., between the "admissibility boundary" specifying the span of the positive- 
definite metrics inside their selected extremal subset and the strength of the 
interaction, respectively. Although the resulting particular curve 7 = 7(e) 
may be specified by a quadruplet of explicit formulae for the segments of its 
boundary, 



7 = 7i ±] = -2± •\/4±2ev / 8-4e 4 + e 6 + 4e 2 -2e 4 , 

Figure [T] offers a better display of all of its relevant features. We should 
emphasize that our considerations are now non-perturbative so that the in- 
dependent variable e G (— v2, \/2) need not stay small. The picture covers 
its full range. We may conclude that the interior of all of the four closed 
loops in Figure [1] represents the prohibited area in which the determinant F 
is negative so that the requirement of the positivity of the metric matrix is 
manifestly violated there. 



14 



4.3 Verification of the ansatz at N = 6 



The A— dependence of all of the six eigenvalues of matrix H^ 6 >(X) may be ex- 
pressed in closed form as well. They all prove real (so that the metric 0^ 6 ^(A) 
exists) for all A G (—1, 1). The metric is obtainable either via its spectral rep- 
resentation [TS1 [TB] or, more easily, from eq. f[2"Uj) . via its computer-assisted 
solution. The resulting matrices 0^(A) form a six-parametric family 



at (1 — A) a 2 (1 — A) 

a s (1 - A) at (1 - A) + a s (1 - A) 

a s (1 - A) a 2 (1 - A) + a 4 (1 - A) 
«4 «5 (1 - A 2 ) + a 5 

a 5 a 4 + a 6 

a 6 a 5 



a 3 (1-A) 
a 2 (l-X) + a 4 (1-A) 
a, (1 - A) + a s (1 - A) (1 - A 2 ) + a 
a 2 {l-X 2 ) + a 4 (1-A 2 ) 
a 3 (1 - A 2 ) + a 5 

a 4 



5 a 

a 6 



(1-A) 



ai 



ot 5 
a 4 + a 6 
a 3 (1 - A 2 ) + a 5 
+ a 5 (l + A) a 2 (1 + A) + a 4 (1 + A) 
'l + X) + a 4 (1 + A) aj (1 + A) + a s (1 + A) 

a 5 (l + A) a 2 (l + X) 

This formula can be read as a confirmation that in our model an optimal 
representation of the general metric will be based on the use of the general 
ansatz 



a 4 

a 3 (1 - A 2 ) + a 5 
a 2 (1 - A 2 ) + a 4 (1 - A 2 ) + a t 
;i + A) + a 5 (l + A)(l-A 2 ) 

a 2 



& N \X) 



N 

£ 

3=1 



a^M) 



(N) 



(A). 



(24) 



Its A— dependence is solely carried by its iV— dimensional sparse-matrix co- 
efficients. The matrix elements (1 ± A) containing the minus-sign always sit 
in the left upper triangle (i.e., above the second diagonal) and vice versa. 
This antisymmetry of the sign of A with respect to the reflection of the ma- 
trix 0( 2X )(A) by its second diagonal simplifies the notation and will hold, 
incidentally, at all the integers K — 1,2,.... 

The individual A— dependent coefficients Mj N \\) are sparse matrices at 
all the dimensions iV = 2K. For illustration, the real and symmetric matrix 



a 3 
a 2 



a 6 
a 5 
a 4 

(1 + A) 
(1 + A) 
(1 + A) 
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(A) may be recalled at K — 4, 

1-A 

1 - A 1-A 2 
1-A (1 - A) (1 - A 2 ) 

1-A (1-A) (1-A 2 ) (1-A 2 ) 2 

1-A 2 (1-A 2 ) 2 

1-A 2 (1 + A)(1-A 2 ) 

1-A 2 

1 + A 

This expression exhibits a seven-diagonal band-matrix structure. The similar 
band-matrix structure will be exhibited by all the matrices M- (A) at any 
N = 2K and at all the subscripts j with the exception of the last one, 
j = N. Precisely this property of the matrix components of the metric (1211) is 
responsible for the possibility of the introduction of a nontrivial fundamental 
length 9. 

5 Extrapolation towards any N = 2K 

Starting from N = 8, matrices Mj N \\) get large and cease to be printable 
easily. Still, the computer-supported symbolic manipulations with these ma- 
trices remain routine and straightforward. Moreover, their structure acquires 
certain features which enable us to extrapolate their low— iV forms to all the 
even dimensions N = IK and test the validity of our extrapolations, with 
much less effort, afterwards. 

One of the most important extrapolation tricks involves the observation 
that all of the matrix elements of metric components (A) form just the 
collection of the following sequence of polynomials 

P = l, Pf ) = l±\, P 2 = l-A 2 , P 3 (±) = (1 ± A) (1 - A 2 ) , 

P 4 = (l-A 2 ) 2 , P 5 (±) = (1 ± A) (1 - A 2 ) 2 , P 6 = (l-A 2 ) 3 , .... (25) 
Their allocation is also not too difficult. 



5.1 Indexing arrays 

Let us consider the mapping Mf°(A) ^ S { k N) between our matrix expan- 
sion coefficients and certain arrays of the same size. This specifies, uniquely, 
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each A— dependent matrix coefficient in series 
At the simplest choice of N = 2 we have 

M{ 2) (A) = 
M 2 (2) (A) = r', <=> .S' 



using an auxiliary array. 








p<+) 









(2) 



(2) 











(26) 
(27) 



At N = 4 this offers a method of an easy coding or reconstruction of the 
four matrices (T23]) . proceeding via the following four indexing arrays Sf at 
j = 1,2,3,4, respectively, 



1 



1 



1 2 
2 1 
1 





1 
1 












The same observations can be formulated at A^ = 6 where the indexing arrays 
form the following sextuplet sf^ at j = 1, 2, . . . , 6, respectively, 



1 



1 1 
1 2 
2 1 
1 1 
1 



1 2 

3 2 

2 3 1 
2 1 

1 





1 
1 2 
2 1 
1 































1 
1 



Once we summarize these N = 4 and A^ = 6 computer-generated results as 
well as their N = 8 and A^ = 10 descendants we reveal the existence of the 
following universal rules. 

• For any polynomial entering any matrix element m\ N \\) the 

L J Jik 

superscripted ± sign must be chosen as + when i > k or as — when 
i < k or as absent when i = k. 
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The indexing symbols as denned by eqs. fl26|) and fl27|) at k — 1, 2 
and N = 2 get generalized to any dimension iV = 2fT. They always 
contain either empty entries or non-negative integer entries "n" . 



• The numerical value of each entry "n" must coincide with the value of 
the subscript n of the related matrix element Pn in the corresponding 
matrix M { k N \\). 

In the light of these rules the complete determination of the functions (A) 
[needed in formula requires just the knowledge of the related indexing 
arrays Sf°. The decoding S^M using the above three rules would enable 
us to reconstruct all the metric matrices Q^ 2K \\) via formulae (12~S1) and (I2"5l) . 
At any .K" the indexing arrays S^ 2 ^ with arbitrary subscript j = 1,2,..., 2K 
can be computed in recurrent manner. The description of the details of such 
a recipe will be provided in the rest of this section. 

5.2 Recurrences for off-central indexing arrays S^ K \ 



3^K 



The inspection of the symbols Sf evaluated by the direct methods at 
the first few integers K = 1,2,... reveals that at any given N = 2K the 
explicit form of the first K — 1 matrices Sj 2 ^ with j = 1, 2, . . . , K — 1 may 

immediately be deduced from their predecessors Sf K ^ . The core of such a 
recurrent construction consists in an enlargement of the dimension followed 
by a symmetric attachment of the two j— plets of units "1" in the empty 
parts of the left upper corner and of the right lower corner. 

The last K matrices S^ C J 1 _j with j = 1,2, ... , K become formed in 

similar manner. Their K predecessors S^^^Jj must be modified by attaching 
j zeros "0" in the right upper corner and in the left lower corner. 

In both these "leftmost-subsequence" and "rightmost-subsequence" sce- 
narios, the results displayed in section 15.11 offer a sufficiently instructive il- 
lustration of the recipe. They also indicate that at the "central" subscript 
j = K the construction of the most complicated missing member of the 
family must be discussed separately. Although it naturally belongs to the 
"leftmost" subsequence, its (2K — 2)— dimensional predecessor (to be denoted 
as proves different horn, the naively expected matrix S^ K . 
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5.3 Recurrences for central indexing arrays S% K ^ 

The sequence of the "middle" or "central" matrices S^ K ^ should be treated 
as exceptional though still generated by a recurrent recipe. Its idea will 
rely on the use of specific predecessor matrices C^ 2K ~ 2 \ Of course, we shall 
proceed in the "leftmost-subsequence" manner enlarging the dimension of C 
and filling K units "1" in the left upper corner and in the right lower corner. 

In order to define the suitable predecessors £( 2K ~ 2 ^ let us start from the 
old "middle" matrix S^5i ^ and apply a specific two-step recipe. Firstly 

we replace each "old" numerical element in S^[ 2 ^ by its successor, i.e., we 
replace "old 0" by "1", "old 1" by "2", etc. In the second step we form 
a left-right reflection of the resulting matrix and arrive at the final form of 
the necessary predecessor £( 2i ^ 2 ) as a result. Thus, at N = 4 we have the 
sequence 

2 



SI 



(2) 



" 1 




" 2 


1 




2 



£(2) 



The recurrent reconstruction of the matrix S^ 2 at iV = 6 will now result 
from adding six units "1" to the auxiliary predecessor matrix in the 
formula 



8P = 



1 2 
2 1 
1 



2 

2 3 
3 2 
2 



2 

3 2 
2 3 
2 



: (6) 



etc. We may conclude that the central matrices S^ K ^ at the respective 
K — 1, 2, 3, 4 (etc) form the sequence 



1 

1 2 
2 1 
1 



1 2 

1 3 2 

2 3 1 
2 1 

1 



1 2 
13 2 

13 4 2 
2 4 3 1 

2 3 1 
2 1 

1 



(etc). The general pattern of their recurrent production is obvious and our 
mathematical construction is complete. 
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6 Discussion 



Our main mathematical result (124"|) represents the general metric matrix as its 
expansion in terms of its sparse-matrix components Mj (A). The latter set 
has unambiguously been determined by certain elementary indexing arrays 
Sf } . Let us now turn attention to some consequences of such a result in the 
context of physics. 

6.1 Fundamental length 9 at large N 

Once we pick up the virtually trivial special case where A — > our present 
toy-model Hamiltonians H^ N '(0) become Hermitian, i.e., tractable as oper- 
ators of an observable quantity also in TiS*>. Then the elements \ip) of this 
friendly Hilbert space may also be re-interpreted as acquiring an entirely 
standard probabilistic interpretation. This leads to a rather exotic arrange- 
ment (published in [18j) where the strict free- motion analogue of our present 
A = Hamiltonian has been assigned the family of non-standard metric 
operators at iV = oo, 

Q(Mostafazadeh) _ j . _ j, . _ ^g) 

Here, symbol X represents the identity operator while V denotes parity and 
k is an arbitrary real constant (cf. eq. Nr. (64) in loc. cit.). The absence of 
any fundamental length has been postulated during the derivation of eq. ff28|) 
in [18]. Thus, in our present Runge-Kutta discretization scenario the latter 
result may simply be interpreted, via identifications I — > M± (0) and V — > 
Mjy (0), as a particular choice of parameters with vanishing 02 = 03 = 
. . . = ajy-i = and nonvanishing a\ 7^ and a at 7^ 0. Thus, in our present 
terminology one has the vanishing elementary length 9 = at k = and its 
unbounded alternative 9 = 00 at n 7^ 0. 

Let us now assume, in contrast, the existence of a finite fundamental 
length in our schematic square-well-type example where the total width 2L 
of the well need not vary but where the dimension iV = 2K itself must always 
grow to infinity in continuous limit. In a preparatory step let us assume that 
N is fixed and that our toy-model metric O^(A) is a (2R + 1)— diagonal 
band matrix. In such a case the natural elementary-length candidate is 
9 = 2Rh = ALR/(N + 1) ~ R. 

In the next step where N starts growing the situation becomes more 
complicated because the elementary-length candidate 9 > should be kept 
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iV— independent at large N. We arrive at the conclusion that whenever we 
try to fix 9 ^ and perform the limiting transition h — > (i.e., A" — > oo), 
we must let the number ~ 2R + 1 of non- vanishing diagonals grow with the 
dimension. This is the reason why all the metrics must be available. 

This observation enhances the relevance of our present solvable model 
where all the necessary constructions were exact. In some purely phenomeno- 
logical applications of the theory with nontrivially nonlocal metric we may 
just search for some approximate results and keep the lattice spacing h fixed. 
Then we are allowed to fix the value of R and to select and construct just the 
particular subset of metrics with 2R + 1 diagonals. Even in such an entirely 
pragmatic setting our present oversimplified example (which admitted the 
changes of both N and R) might still prove useful as a methodical guide. 

Alternatively, we might pick up a nontrivial R > and keep this integer 
(i.e., the number of diagonals in 9^) fixed even during the limiting transi- 
tion A" — > oo (yielding, formally, 9 = of course). Then we still obtain the 
sequence of metrics which may converge, typically, to a nontrivial operator 
@(oo) represented by a generalized, momentum-dependent kernel entering an 
appropriate operator generalization of our present, smooth normalization in- 
tegral etc. Naturally, this is a promising but mathematically difficult 
possibility which we couldn't have addressed here. 

6.2 The reality of the energies at large N 

From our most elementary toy Hamiltonian H^ 2 '(X) we deduced the energies 
most easily, E = E { ± ] = 2 ± y/l - A 2 . They remain real in the closed interval 
of A G (—1,1). The purely numerical analysis of a few further members 
H^ N '(X) of the family reveals that the related energy spectra remain real in 
the same interval of the couplings A G (—1,1). 

Empirically, the validity of this observation is illustrated here by Figures [2] 
and[3j Sometimes, it may prove useful to re-parametrize A = cosip G (—1, 1) 
with ip G (0, 7i"), therefore (cf., e.g., paragraph 13.31 above). 

One should add that several other features of bound states (e.g., the eval- 
uation of matrix elements of some operators of observables) remain transpar- 
ent and well illustrated by our discrete short-range model H^ N \\) and, in 
particular, by its most elementary special cases (fTT]) - (|T3l . 
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Figure 2: Spectrum of H^(X). 



6.3 The role of VT— and CVT— symmetry 

On the background of the existence and undeniable physical appeal of the 
so called VT— symmetric differential-equation models as reviewed by C. 
Bender [10J one reveals that also all the elements of our sequence of the 
Runge-Kutta discretized models can be incorporated in the same context. 
Indeed, our Hamiltonian matrices ( ITTi) . ( fl2l) . ( fl3l) etc may be identified as 
VT — symmetric, provided only that we treat the operator T as mediating 
transposition and that the parity V is represented by the matrix with units 
along its second diagonal, 



The specific merits of our discrete model involve, furthermore, the simplic- 
ity of mathematical analysis since we avoided the Rayleigh-Schrodinger-type 
perturbation expansions reported, in refs. J2TJ [18], as fairly difficult and 
complicated. In this context the Runge-Kutta discretization opened a way 
towards our straightforward and efficient application of linear-algebraic tech- 
niques. 

The post-multiplication or pre-multiplication of any matrix H by the par- 
ity matrix V mediates the left-right or up-down reflection, respectively. We 
immediately see that our toy Hamiltonians are not only VT — symmetric but 
also, in the standard terminology of linear algebra, V— pseudo-Hermit ian 



*Pi,n — ^2,AT-l 



1 . 
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Figure 3: Spectrum of H^(X). 

[T5] and G— quasi- Hermitian [19]. In this spirit we may recollect ref. pU] 
and try to factorize the metric, = CV. Without the usual additional 
constraint C 2 = I this merely defines certain additional, "protocharge" oper- 
ators C. Vice versa, the incorporation of the constraint C 2 = I converts the 
protocharge factors into the popular "charges" jTUl [31]. In principle, their 
existence imposes a fairly severe constraint upon our freedom in the choice 
of the parameters in the metric [32]. At the same time, these constraints 
may still be expected to leave some residual freedom in the domain T> of 
admissible parameters aj [TH [33] . 

7 Summary and outlook 

In the mathematical part of our paper we presented an exact construction of 
the most general metric operator 9^'(A) (i.e., of the most general positive- 
definite inner product) for one-parametric non-Hermitian matrix Hamiltoni- 
ans H( N \\) introduced in paper [TTj . The construction proceeded in two 
steps. Firstly, the brute-force use of computer enabled us to generate an 
exhaustive list of all the admissible matrices ®( N >(\) at the first few (even) 
matrix dimensions N = 2K. In the second step we revealed a clear pattern 
of the dependence of these low-dimensional matrices on their dimension N 
and on the coupling constant A. This enabled us to arrange each Q^ N \X) as a 
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linear superposition of its N elementary sparse-matrix components M- (A) 
with band-matrix structure. In the third step an ansatz for the latter com- 
ponents has been found and its validity has been confirmed by extrapolation 
and its subsequent facilitated verification. In the final, fourth step the deter- 
mination of the matrix elements in Mj (A) was reduced to their indexing via 

arrays Sj with integer or empty entries defined via an elementary recurrent 
recipe. 

Beyond the horizon given by our particular illustrative example we paid 
our main attention to the rather serious problem of the constructive approach 
to the models with fundamental length 9 and, in particular, to the practical 
feasibility of the necessary construction of the related, "fine-tuned" metric 
operators 0. We tried to overcome the well known difficulties encountered, 
in the literature, during perturbation constructions of the metrics. We found 
a way how to get rid of at least some of the current methodical obstructions 
resulting from the immanent weakness of perturbation techniques. In this 
context, the so called Runge-Kutta approximation techniques were found 
very productive and strongly recommendable. 

In the parallel physics-motivated discussion of the relevance of our results 
in the abstract quantum theory as well as in its various applications we 
remind the readers, first of all, that our toy Hamiltonians H^ N '(\) with iV = 
oo and A 7^ did already serve as a guide during the recent discussion and 
clarification of some conceptual problems concerning the quantum scattering 
by non-Hermitian point interactions [11]. In the present continuation of their 
study we succeeded in reconfirming the relevance of the similar schematic 
interaction models also in the theory of bound states. 

Our main attention has been paid to conceptual questions again. We 
opposed, e.g., the frequently postulated absence of fundamental length in 
VT— symmetric and non-Hermitian models. The resolution of some con- 
crete technical problems has been found. In particular we revealed that the 
extreme simplicity of our model opens an interesting nonperturbative way 
towards an innovative and fully constructive understanding of the emergence 
of an elementary length in the quantum system in question. We were also 
able to add a few new ideas to the lasting discussions concerning the interpre- 
tation of the presence of a fixed scale 6 > in quantum theory. We proposed 
that in the language using the concept of metric operators many theoret- 
ical considerations may remain feasible even when a nontrivial constant 9 
is introduced. Last but not least we offered a few arguments supporting 
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the possibility of using alternatives to the popular and widespread strategy 
which connects the "coordinate-smearing" quantity 9 > 0, indirectly and 
exclusively, to a hypothetical non-commutativity of classical coordinates. 
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Appendix I: Sparse matrices M^ N \X) at A = 



The formal structure of the complete sets of metrics 0^(A) collected at 
the smallest dimensions N < Q gets fully transparent in the vanishing- 
potential limit A — > when all the Hamiltonians become Hermitian in the 
iV— dimensional Hilbert spaces TC^ F \ The textbook choice of the metric looks 
unique because people tacitly assume that there exists no nontrivial funda- 
mental length in the theory (cf., e.g., ref. [18J). In our present notation such 
an option coincides with the special case where 9 = 0, a,\ > (i.e., say, 
oti = 1) and a2 = (X3 = ■ ■ ■ = = 0. 

Whenever we intend to build the theory where the choice of 9 > sets a 
nontrivial length scale, Hamiltonian H need not be non-Hermitian in 7i^ F \ 
In our toy model at A = 0, in particular, we may still define a non-Dirac 
metric using formula (1241) with some nonvanishing values of parameters «2 
and a?3 etc. It is only necessary to guarantee that the metric matrix itself 
remains positive (for illustration, recollect inequality (|19p which specifies the 
full allowed range of parameters at N — 2). 

Although our models H^ N \X) admit a nontrivial fundamental length 
9 > even in their square- well limit A = of paper [21], possible physics 
represented by such an extreme example looks rather artificial, Still, its me- 
thodical merits are remarkable. Firstly, the coefficient matrices M = M,- (0) 
become solely filled by the matrix elements or 1. Secondly, their knowledge 
may prove useful for coding the indexing arrays in computer-assisted manip- 
ulations. Thirdly, the simplicity of the model implies that the j—th member 
of the sequence can be defined by the following closed formula, 

(Mf^ , fc (°) = 1 iff i-k = m, N + l-i-k = n, (29) 

m = j- l,j-3,...,l-j, n = N -j,N -j - 2, . . . , j - iV 

and the verification of validity of this formula by its direct insertion in eq. (j20p 
is very quick. Fourthly, the existence of this and related formulae may prove 
useful for perturbation constructions in weakly non-Hermitian dynamical 
regime where |A| 1. 
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Appendix II: Our toy Hamiltonians at large N 



At a sufficiently large N our particular one-parametric Hamiltonians H( N '(X) 
may be interpreted as discrete versions of a differential operator with a point 
interaction localized in the origin. For a deeper understanding of such a 
correspondence let us abbreviate 2 — h 2 E = cos e as usual [21]. We may then 
treat e G (0, ir) as a new energy variable and visualize the wave functions 
i/j(x) with x ^ as satisfying a free-motion equation complemented by the 
respective left and right initial conditions r/>(— L) = = 0- At a fixed L 
and in the N = 2K 3> 1 approximation these free-motion-like solutions must 
be further restricted by the pair of A— dependent constraints near the origin, 

(1 + X)Mi) -2cose^(a; ) + M-i) = 0, (30) 

ip{x 2 ) - 2coseip{x 1 ) + {1 - X) ip{x ) =0. (31) 

In the limit h — > we may expect the emergence of a discontinuity in i/j(x) 
at x = 0. Even at all the finite N ~ 1/h ^> 1 the wave functions remain 
well represented by their respective one-sided Taylor series near x = so 
that eqs. (1301 and fl3Tl) may be interpreted as a matching condition. Once 
we return to the original energy variable h 2 E = 2 — 2 cose = F and insert 
the truncated expansions 

V(x_ x ) = Vfc(0) - |^V^(0) + 0(h 2 ) , ^ ) = 1^(0) - \hij' L (0) + 0(h 2 ) , 

Tp(xi) = Mo) + \h^'M + o(h 2 ) , *p{x 2 ) = mo) + \hMo) + o(h 2 ) 

in eqs. (1311 and fl3"T|) . a straightforward algebra leads to the following ele- 
mentary condition 

hf-(l + X) F+l \f M0) \ = ( 1 + A F-1W^(0)\ 

2 ^-(f+i) i-a;^l(°)/ v^- 1 i-aA^(°); 

(32) 

which matches the wave functions and their derivatives in the origin. 

In the domain of sufficiently small h > the latter relation is equivalent 
to the original constraints (l30j) and (!3T|) . We may conclude that at all the 
nonvanishing small h > our conditions (|32j) leave our interaction in the 
origin translucent and manifestly energy-dependent. 
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Various special cases of our N ^> 1 bound-state model may be studied 
noticing, for example, that the energy- dependence disappears in the low- 
excitation regime where the quantity F = h 2 E remains negligible. At a 
generic energy F > the above set of h > solutions must be comple- 
mented by the two additional, anomalous bound states emerging at the two 
exceptional energies F = F± = 1 ± yl — A 2 which make the coefficient matrix 
singular. At these energies, both the values of ?Pr } l(0) are, in general, non- 
vanishing and firmly determined by our choice of the two derivatives ip' L {—L) 
and ip' R (L) at x = =fL, respectively. Thus, both the left and right branches 
of our two exceptional bound states are obtained by the same matching in 
the origin as above. Their specific feature is that at both our exceptional 
energies F± the two lines of eq. (I3"2|) degenerate, at all the sufficiently small 
h ~ 0, to the sing le constraint ip R (0) vTTA ± ip L (0) y/1 - X = 0. Thus, 
manifestly asymmetric wave functions are obtained. 

In the continuum limit iV — > oo our sequence of the matrix Hamiltonians 
H^ N '(X) can be reinterpreted as a series of dynamical models which converge 
to a specific differential equation with a point interaction potential in the 
origin. It is easily seen that at a generic energy E the h — > limit of eq. (J32]) 
leads to the vanishing F = 0(h 2 ) so that the above-mentioned "exceptional" 
solutions disappear from the spectrum. At any A ^ 0, only the elementary 
opaque-wall constraint i/}r(0) = t/>l(0) = survives and leads, say, to the 
two independent series of bound-state solutions which live solely on the left 
or right half-interval of x, respectively. 
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